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Summary 


Experiments have been carried out to test the accuracy of laser-Doppler 
instrumentation for measurement of Reynolds stresses in turbulent boundary 

” layers in supersonic flow. Two facilities were used to study flow at constant 
pressure. In one facility, data were obtained on a flat plate at M, = 0.1, 


7 with Re, up to 8,000. In the other, data were obtained on an adiabatic nozzle 


6 


wall at M, » 0.6, 0.8, 1.0, 1.3, and 2.2, with Re, = 23,000 and 40,000. The 


cs) 
present Part I of this repert is limited to a description of the mean flow as 
observed using Pitot-tube, Preston-tube, and floating-element instrumentation. 
Emphasis is on the use of similarity laws with Van [cieet scaling and on the 
interence of the shearing-stress profile and the normal velocity component from 
the equations of mean motion. The experimental data are tabulated. 

Part II of this report, published separately, is a description of the 


mean flow and Reynolds-stress field as observed in the same f'ows using ¢ 


single-particle laser-Doppler velocimeter. 
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I. Intredsction 
The turbulent boundary layer at constant pressure has been studied 

experimentally and theoretica!ly for many years. Collected experimental 
mean-velocity data for low-speed flow have been carefully reviewed by 

Coles (1962, Appendix A), who recommends taking the measurements by 

Wieghardt (1943) as the best available standard. A catalog and a comparable 
review of mean-velocity data for high-speed flow (including flow with pressure 
gradient) are presently being prepared for AGARD by Fernholz and Finley (1977). 
One method being used by Fernholz (1976) to organize the information in this 
AGARD catalog is recasting of the compressible-flow results into a form 
appropriate for incompressible flow. For this purpose, the wall-wake model 

for the mean-velocity profile (a model which has been thoroughly exercised 

for incompressible flow by Coles (1968)) and the mixing-length scaling proposed 
for compressible flaw by Van Driest (1951) appear to be quite useful. 

Measurements of Reynolds stresses in high-speed turbulent boundary layers 

are rare. For incompressible flow, the turbulent shearing stress can be mea- 
sured directly, or it can be calculated from the distribution of mean velocity 
with the aid of well-established similarity laws. Good agreement between mea- 
sured and calculated values, as in the case of the hot-wire measurements by 
Klebanoff (1954), helps to establish confidence in the extension of hot-wire 
methods to more complicated flows. Recent measurements by Johnson and Rose 
' (1973), Yanta and Lee (1974), and Abbiss (1976) have attempted to extend this 
process to the use of laser-Doppler instrumentation in supersonic flow at Mach 
numbers in the range 1.5 to 3.0. However, a serious anomaly appears in the 
case of the turbulent shearing stress, defined as =pu‘v’. The maximum value 


occurs much further from the wall than is reasonable for flow at constant 
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pressure. The anomaly has been discussed by Sandborn (1974), who supports 
the conjecture by some of the authors cited that density fluctuations may 
contribute substantially to the turbulent stresses near the wall. This 
conjecture is in direct opposition to the conclusion by Morkovin (1961) that 
effects of density fluctuations should be small compared to effects of vari- 
ations in mean density for Mach numbers uf to 4 or 5. 

The purpose of the present experiments is to obtain redundant data 
over a substantial range of Mach numbers (M, = 0.1 to 2.2), in an effort 
to resolve the anomaly in turbulent shearing stress. Essentially the same 
range of Mach numbers has also been studied by Winter and Gaudet (1970), who 
used a Pitot tube to determine mean velocity and a floating-element balance 
to determine surface friction. However, no measurements were made of turbulent 
stresses, 

The present measurements have sufficient redundancy to permit a realistic 
assessment of their accuracy. The primary instrumentation is a Pitot tube which 
traverses the boundary laysr. In addition, surface-friction measurements are 
made using both a float‘ng-element balance and a Preston tube. The mean-flow 
scaling suggested b:; Van Driest is applied to the data, to test the adequacy 
of a single similarity formulation for both compressible and incompressibie 
flow, and the shearing-stress distribution is calculated as part of the analysis. 

The results are discussed in Part I of this work. An Appendix contains 
a complete record of the experimental data in tabular and graphic form. 

A laser-Doppler velocimeter has also been used to measure mean velocity 
and three components of Reynolds stress in the same flow. The results of the 
LDV studies will be reported by P. E. Dimotakis, D. J. Collins, and D. B. Lang 


in Part II of this work. 
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Il. Flow Facilities 

A. High-Speed Flow 

Measurements were made in the ceiling boundary layer of the 20-inch 
wind tunnel at the Jet Propulsion Laboratory, at nominal free-stream Mach 


numbers M, of 0.6, 0.8, 1.0, 1.3, and 2.2, at nominal Reynolds numbers Reg 


of 23,000 and 40,000. The JPL facility is a continuously operating, variahle- 
density tunnel, with a test section 45.7 cm wide by 50.8 cm high. The top 
and bottom walls of the tunnel diverge slightly to compensate for boundary- 
layer growth. 

For the present experiments, the region of uniform flow in the test 
section was extended approximately 150 cm beyond the end of the flexible 
nozzle, or 60 cm beyond the center of the schlieren windows, by installation 
of a pair of instrumented flat plates on the floor and ceiling of the tunnel. 
Particular care was taken to obtain a smooth junction between the plates and 
the nozzle wall. 

From the experience, for example, of Liepmann and Ashkenas (1947), it 
has lony. been known that the experimental treatment of the downstream boundary 
condition is important at transonic speeds. The unsteady behavior often ob- 
served in transonic shock-wave boundary-layer interactions may be partly a 
consequence of unsteady flows generated in the diffuser. For tne present 
experiments, stable flow at high subsonic free-stream Mach numbers was achieved 
by introducing a variable-thickness double-diamond airfoil choke in the dittuser.” 
The choke was oriented vertically in the diffuser, normal to the test plate, with 


the leading edg» of the choke located 70 cm downstream from the balance station. 


* 
The advice of H. Ashkenas during this development is much appreciated. 
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The point of waximum thickness was located 15 cm further downstream. Flow 
past the choke was religved by expanding the diffuser doors by 7.3° on each 
side in order to maintain constant area at the minimm thickness setting. 
High-speed schlieren movies showed that this arrangement eliminated the 
upstream-running waves observed in previous experiments and provided a 
quiet environment in which to perform transonic boundary-layer experiments. 

Local static pressure was measured at 82 static-pressure orifices 
located throughout the test section and diffuser. The measurements used 
the JPL miltiport measuring system, which simultaneously recorded the stag- 
nation temperature and pressure, the free-stream static pressure, and the 
pressure from two 0-15 psia Statham pressure transducers, each of which 
sequentially sampled 50 orifices. 

Typical free-stream Mach~number distributions for the present experi- 
ments are shown in Fig. 1. There is no substantial pressure gradient over 
a distance of about 140 cm upstream and about 40 cm downstream from the 
balance station. 

B. lLow-Speed Flow 

Additional measurements were made in the Merrill wind tunnel of the 
Graduate Aeronautical Laboratories at the California Institute of Technology, 
#.¢ a free-stream ve ‘ocity of 37 m/sec. This tunnel is a continuously operating 
closed-return facility with the downstream end of the test section vented to 
ambient pressure. The test section is 115 cm wide by 82 cm high and has diverging 
walls to account for boundary-layer growth. 

The test plate for these experiments was made from 1.9-cm thick plywood, 
surfaced on both sides with l-mm thick formica to provide a smooth finish. 


The leading edge was elliptical, with a transition strip located immediately 
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downstream from the elliptical section. The horizontal plate spanned the 
test section and extended 244 cm downstrcam from the beginning of the test 
section. The plate was supported from the ceiling of the tunnel. and all 


measurements were tarde on the lower surface. 


Twenty static pressure taps were provided on the surface cf the plate. 

A Scanivalve was used to select the pressure to be read by a Barocel digital 
manometer. The resulting free-stream Mich-number distribution, shown in 
Fig. 2, indicates that pressure-grudienct effects should be small. 

TIT. Pitot-Pressure Data 

A. Instrumentation 

For tne experiments at JPL, a Pitot-pressure probe could be introduced 
into the boundary layer through the ceiling of the tunnel at any one of the 
five axial stations listed in Table 1. The origin for the x-coordinate is 
the center of the floating-element friction balance, 11.3 cm downstrear from 
the junction between the nozzle wall and the test plate. During the probe 
measurements, the balance vas replaced by a blank port which was instrumented 
with static-pressure taps. 

The Pitot-pres ~re probe was constructed from stainless steel hypodermic 
tubing. The probe tip was formed by flattening 0.127-cm diameter tubing to an 
oval measuring 0.0127 cm inside (in the direction normal to the plate), with the 
lip thickness honed to 0.003 cm. The center of the support stem was 5.08 cm 
downstream from the probe tip. The probe position, the Pitot pressure, the 
tunnel stagnation temperature, and the tunnel stagnation and static pressures 
were recorded by the data system. 

For the experiments at CIT, two techniques were employed. Within the 


first 100 cm from the leading edge of the plate, Pitot measurements were made 
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using a seven-tube rake. Further downstream, Pitot measurements were made 

by traversing a small probe through the boundary layer, as in the high-speed 
experiments. The probe tip was flattened to an oval measuring 0.0203 cm inside 
(in the direction normal to the plate), with a lip thickness of 0.020 cm. 
Boundary-layer measurements were made at the stations listed in Table 2. 

B. Data Reduction 

For each Pitot-pressure profile, a change of slope in the pressure data 
was used to uefine the point of contact of the probe with the wall. No dis- 
placement correction was made. The free-st~cam static pressure for each profile 
was taken as the average static pressure in the te section in the vicinity of 
the probe. The flow properties at the edge of the boundary layer were then 
computed using the average Pitot pressure well outside the boundary layer. 
Assuming constant static pressure, the local Mach number was computed either 
directly or through the normal shock relations, as appropriate. 

The local stagnation temperature in the boundary layer was not measured. 
However, for the JPL experiments, the temperature measured by a thermocouple 
embedded in the surface-friction balance structure indicated that the flow was 
essentially adiabatic. Hence the temperature may be estimated from a variant 


: , * 
of the adiabatic Crocco relation, 


*Equation (1) is often used in the reduction of experimental data, despite the 
fact that it does not conserve energy in adiabatic flow. The present data 
analysis assumes that p and M are measured exactly. Hence so are Ti/t bs 

1+ (y-1)92/2, u/CyRT) ® = M, and Pt: = yp. If use of Eq. (1) introduces 

a local relative error of ¢ in To» the relative errors in T, u, 0, and pu 


are ¢, e/2, -€, and -c/2, respectively. 
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6 4 2 1 + (X54) u? 
T vy-1l 2 u r 
ies: a = ltr ( ) M | = (2 ) | = pete ) 


where the recovery . tor r is defined by 
a - T, 
| a 2) 
T -T 7 <2) 
° e 
e 


and is assigned the constant value ~ +: 0.885. 

C. Results for the Mean Flow 

Typical mean-velocity profiles measured at the balance station 
(JPL-4) for nominal Reynolds numbers Reg of 23,0C0 and 40,000 arc presented 
in Figs. 3 and 4. Values for viscosity are obtained from the Sutherland 


viscosity law, 


$ 3/2 
E o(H= eye)” (3) 
vy T+S/\T , : 


where T_ = 291.75 °K, S = 110 °K, and uy, = 1.827 x 1974 


gm/cm-sec. One 
profile at Re, = 8000 f-0m the low-speed experiments (CIT-9) is a1.0 
included in the figures for comzirison. A complete data tabulation 


appears in the Appendix. 


{ntegral thicknesses for the boundary layer are computed from 


.) 
* ou 
os L- S—]}a4y, (4) 
(0) Peve 
and 
L) 
o- | fu (1-2 ey (5) 
u 
lo ee e 
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The boundary-layer form parameter H is defined as 


= 
« 
O|7, 


é (6) 


For two-dimensional mean flow, the surface friction can be obtained 


from von Kdrmin's womentum-integral equation, 


dp 2) 9 1 aP 
Cc, 2 $2.2 (2+n- 0,7) 2p pe (7) 
e 


The accuracy of Eq. (7) is expected to be low, ,rimarily because of dif- 
ficulty in differentiating experimental data for @{x) and ui Go (see Table 
A3 of the Appendix). For the present measurements, the second tm in 

Eq. (7) is at most 3 percent of the first term, and is uncertain by a 


comparable amount. Hence this term has been ciscarded. Values for 


Ce = 2 dA/dx are listed in Table 3, which compares values obtained for 


C. by this and several other methods. 


D. Van Driest Scaling 


The compressibility transformation proposed by Van Driest (1951) 


uses the mixing-length expressioa 
2 
tTeT - ot? (22) (8) 


. together with Prandtl'’s hypothesis 


Lenny, (9) 
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to obtain 


1/2 du z os 


dy ny (10) 


1 
The appearance of the combination (p Pau) Sugyests that the velocity u 


should be replaced by an effective velocity u™ def ined by 


u 
sous wae pace 
ui = (=) dui, (ity 


Q 


where o* is a constant reference density included for dimensional reasons. 
Integration of the mixing-length equation (8) then gives 


1/2 
* if y. 
eS (4 ln (= +constant , (12) 


where y* is a constant reference length also included for dir- nsional 
reasons, 

Equation (12) is typical of mixing-length formulas in that it is at 
best er unclear description of a z.nall fragment of the mean-velocity 
profile. The choice for o* and y* and the value of the constant in 
Eq. (12) are customarily resolved by emphasizing quantities evaluated at 
the wall. For example, the definition (11) is readily integrated in closed 
form for the energy integral (1). The result is the Van Driest scaling for 


velocity in the case of adiabatic flaw, 


i. 4 
m (=) = = cin’ (a = ) ; (13) 
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where m, defined by 


m ae ree (lym , (14) 


obviously cannot exceed unity. 


The form of Eqs. (12) and (13) suggests, but does not require, 


= * = wos 
choosing o o, and y vies where 
(3) 
ule) ’ (15) 
7 o, 


is the friction velocity. The choice y* = v/a, s in particular, is 


necessary if the functional dependence of u on y in Eq. (12) is to hold 


at the wall. Such reasoning, however, is not part of the mixing-length 
argument, which applies only outside the sublayer. Given these choices, 


then in a usual notation Eq. (12) becomes 


ae = : ln y* +e , (16) 
whe re 
+ u* + Yr 
us = ee Galas (17) 
7 w 
and 
u* -l u 
ns = sin ( v.) . (18) 
e e- 
GE b&b 
ORIGINAL BA 
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The choice for o*, u,> and y* is important because it controls the 
dependence of x and c on M, and y. What is wanted is the particular choice 
which minimizes this dependence. There is substantial evidence, for example, 
ag papers by Fenter and Stalmach (1957), Rotta (1960), Moore and Harkness 
(1964), Maise and McDonald (1908), Michel, Quemart, and Elena (19609), Danberg 
(1971), Squire (1971), and Fernholz (1976), that use of wall quantities as in 
Eqs. (16)-(18) is very nearly optimum from this point of view, at least for 
adiabatic flow at constant pressure at Mach numbers up to 5, 

Most of these authors have also gone beyond the mixing-length argument 
to consider a more general fit to a defect law or to a combined wall-wake 


formulation of the mean profile, in the manner adopted by Coles (1968) for 


low-epeed flow; i.e., a fit to 


a? Pang Fee tae, (iy) 
% xn 
where 
aly 
Y 75° (20) 


Such a fit has been carried out for the present measurements, with quite 
satisfactory results. The constants x and c are given their incompressible 
values, x» = 0.41 and c = 5.0. The parameters uss fl, and 6 are then de- 
termined by a two-parameter least-squares fit of the experimental data to 
Eq. (19), taking as a third condition the constraint imposed by the local 


friction law, 
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As proposed by Coles (1968), data near the wall and near the free stream 
are excluded. For the JPL experiments, data are retained for y 2 200 
and y/8 © 0.95. For the CIT experiments, data are retained for y* 2 8&0 
and y/6 < 0.95. Typical examples of the resulting fit are shown in Fig. 5. 
The values obtained for § are indicated by tick marks in Figs. 3 and 4; 
they correspo:ud to values for u/u, of 0.996 to 0.998, 

The quality of Van Driest scaling, when universal constant values 
are assumed for x and c, can be tested in different ways. One test is to 


compare values inferred for the local friction coefficient 


o, uw : 
e e 


with values obtained by other means. Table 3 makes this comparison. If 

the floating-element data are taken as a standard, the conclusion for the 
presert experiments is that the profile fit gives values for Ce which are 
slightly high. The discrepancy is small at subsonic speeds, but increases 
to about 6 percent at M, = 2.2. 

A second test is to compare values obtained for the profile parameter 

il with corresponding values for low-speed flow, as defined by the low-speed 
data of Wieghardt (1943). This comparison is made in Fig. 6." The main 


conclusion is that there is very little effect of compressibility on the 


shape of the mean-velocity profile in Van Driest coordinates, ..t least for 


* 
The particular choice of CrRe, for the abscissa in Fig. 6 (Coles 1962) is 


8 
not important for these data, although it might become importanc if the 


figure included data at higher Mach numbers and lower Reynolds numbers. 
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Mach numbers up to 2.2. This conclusion is supported by the inclusion 
of a few representative points from the work by Winter and Gaudet (1970). 
‘here may be a tendency for fl to decrease slightly at large Reynolds 
numbers, as noted for low-speed flow by Coles (1962, Appendix A). 


E. Inferred Data for v/u and t/t, 


The distributions of normal velocity and shearing stress through 
the boundary layer are of central interest in these experiments pvecause 
of the direct comparison with LDV measurements to be made in Part II of 


chis report. From the equaticn of continuity, 


y 
3 
wo =| os. (23) 


and from the equation for conservation of momentum in flow at constant 


pressure, 
y ‘ 
pu 
Tt + aera 
y ™ pay f a dy. (24) | 


It is desired to evaluate the integrals in Eqs. (23) and (24) for 
the Van Driest description of the mean velocity profile with similarity, 


Eq. (19). A useful first step is a change of variable. Put 


ux 
mou, (25) 
e 
so that Eq. (28) becomes 
u 
-~« sinU. (26) 
e 
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The corresponding form of Eq. (1) is 


i.) T 
2 
= = = cos U, (27) 
e 


These may be substituted in Eqs. (23) and (24) to obtain 


T 
e d& 
Ove w OMe T. dx Pp (28) ‘ 
w 
and 
T 
2 e db u 
: Ww Mee ty x (20 iz u, P) . (29) 


where the quantities denoted by P and Q are the definite integrals 


y 
2 
ay. { Loan) oY ay: (30) 
md 


dx cos7U ox 
and 
y 
d& 1 sinU oU 
tg. .t, [sit Bay ° (31) 
dx me hy cos2U ox 


Note that dé/dx is a phantom factor in these expressions. If Eqs. (28) 
and (29) are evaluated at the edge of the boundary layer, where Tt = 0, 


O7=p,uu, andv =v, the result is 
e e e 


v T * 

e e d6§ _ a6 

ey dx Fe ~ dx ‘ (32) 
e w 


AEOC-TR-78-21 


and 


T T 
w e df dA 
2 T dx i Page are (33) 
u w 
ee 
It follows that 
T * 
dé w dé 
dx *e "T dx ’ (34> 
e 
and that 
T 
d§ > 1 wd * 
dx a 2 Ty, dx co" + 6) . (35) 


Given U(x,y), the most convenient form for calculation is probably 


the normalized form 


v eP 
Caer aE See (36) 


and 
ee es (37) 
w (7% - Fe) 


To undo the normalization in Eq. (37), a value must be specified for 
1, /0.u, i.e., for Ce: The derivative d&/dx may then be calculated from 
Eq. (33) and inserted in Eq. (32) to obtain a value for v./u,- This value 


can be used in turn to undo the normalization in Eq. (36), with the result 
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‘ (38) 


The analysis so far involves only the formalism of Van Driest 
scaling, inasmuch as the function U(x,y) has not been specified. For 
purposes of curve fitting, this function is defined by Eq. (19) outside 
the sublayer. Other authors, notably Maise and McDonald (1968) have also 
made calculations equivalent to using Eq. (19) in Eq. (29) to obtain the 
distribution of t/t. However, for accurate evaluation of the integrals 
P and Q near the wall, both U and U/3X ne2d to be more accurately defined 


in the sublayer. We therefore revise Eq. (19) to read 


u 
en. + qn 2 
rf U = f(y ) +2 ; sin’Y , (39) 


c 
r) 
Bi~w 


and we describe the flow near the wall by ar. ‘mplicit formula for f(y") 


proposed by Spalding (1961) and independently by Kleinstein (1967), 
+ onc | mf (nf)? Gey? 
y =f€+te [e - 1 - (xf) - a og ] : (40) 


This formula has the proper behavior near the wall, where f = y* + oy*)*, 
and also outside the sublayer, because for (xf) >> 1 Eq. (40) reduces to 
Eq. (16). It is Eq. (40) which is plotted in the sublayer region in Fig. 5. 


For Tl, u, and m constant, differentiation of Eq. (39) gives 


1 Ye au _ [ + + 4 | 
mu, ax LY bi ee a + 


du, I 


Ineo 1 48 
ae 25 ¥ sin(2¥) ¢ ; (41) 


tw 
te 


t 
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where, from Eq. (40), 


2 
st wl tne ** [eré ok SKE). = (ue) ] (42) 
€ (y ) 


The derivative du, /dx in Eq. (41) can be eliminated by noting from Eq. (21) 


that 
du 
1 d8 +, 1 i 
6 dx Dies KU ) up, dx (43) 
Consequently, the integrals P and Q may be written 
1 Yy (x4 Pe + ¢ 
ae | (espe) y £949 ny sin (2y)| ay , (44) 
cos ee a (xu, + 1) 
and 
y 
u + +27 4 
ge ea econ +2 M1 ¥ sin (2¥)| dy . (45) 
eYo \cos"U (wa, + 1) | 


The integrals P and Q are readily determined for a given profile 
once the parameters uss Il, and 6 are specified. For convenience of tau- 
lation, we use experimental values for y and we determine y" from (17), 

Y from (20), U from (25), u* from (39), and £ and £’ by interpolation in 
(40) and (42), respectively. Thus the measured data influence the cal- 
culations only indirectly, through the fit which determines us l, and 6. 

Figures 7 and 8 show typical distributions for v/u and tht, cal- 

culated by this method, using profile parameters taken from the fit de- 


scribed in Section III-D. Frem these figures it is clear that both 
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quantities scale with outer variables, and that in normalized form they 
are relatively insensitive to changes in Mach number and Reynolds number. 
Several authors, including Meier and Rotta (1970), Bushnell and 
Morris (1971), Horstman and Owen (1972), and Sturek (1973), have used a 
different and less structured method for calculating tit, in flow at 
constant pressure. Instead of wall-wake similaricy with Van Driest scaling, 
the basic assumption is that u/u, and p/p, are functions only of y/& or y/@. 
In the former case, the problem of defining § must be faced at the outset. 
In the latter case, this problem can be postponed. For generality, we take 


the independent variable as y/L. Then Eq. (28) is replaced by 


/L 
ov 
yeu de(z “ee fe at), (46) 
dx (L ou ou L 
io ef & 
and Eq. (29) is replaced by 
y/L y/L 

ser - pu? He (et fH at- mg (47) 

w PoYe dx \u ou L 2° Lj’ 

e ee O46 


When the integrals extend to the free stream, these become 


8* dL 
yo wT ax ‘ (48) 


and 


(49) 
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The last two equations are clearly not compatible if L is che 
same for both. To satisfy the condition ve fu, = dé*/dx, it is necessary 
to take = 8°. To satisfy the condition rio,” = dO/dx, it is neces- 
sury to take L = @. In neither case ts L = § a satisfactory choice. 

To illustrate the problem, soma typical results ac ording to 
these equations, with L = 6, are cumpared in Figs. 9 and 10 to earlier 
results based on Eqs. (44) and (45). Experimental points now have a direct 
influence on the calculation, because they define the functions to be 
integrated. The distributions in the figures are therefore properly 
rounded in the vicinity of the boundary-layer edge, avoiding the corner 
which is present in the earlier results. There is a slight problem witt. 
Eqs. (46) and (47) at small values of y/@, where the experimental values 
of M/M,, v/u,s and o/o, are all larger (perhaps because of probe inter- 
ference) than the values associated with the profile fit. The integrals 
thus become permanently biased during the passage through small values of 
y. However, the main source of the discrepancy in v/u outside the bovndary 
layer is the fact that Eq (48) re uires v,/u, = (8*/@)d@/dx, rather than 
the correct value dé*/dx. When the difference @di!/dx is estimated independ- 
ently, using the tabulated material of the Appendix, the discrepancies in 
Fig. 9 are quite well accounted for. Because of these discrepancies, and 
because the wail-wake fit provides an unambiguous definition for 6, we 
consider the calculation based on Eqs. (44) and (45) to be superior. 

IV. Surface-Friction Data 
A. Floating-Element Balance 
A floating-element balance used by Coles (1953) was recommissioned 


* 
for use in the present experiments. The only important design change was 


* 
The expert assistance of George Tennant in preparing the talance is gratefully 


acknowledged. 
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in the method of achieving the null position for the element. The balance, 
shown in Figs. ll and 12, is a sealed unit mounted in a 23.5 cm-diameter 
port flush with the ceiling plate in the 20-inch wind tunnel. The original 
element occupies a 10-cm diameter circle which is iocated 5.08 cm upstream 
of the port centerline, as indicated in Fig. 12. 

The balance is a null device with the floating element supported 
by a four-flexure linkage. The total force on the element is inferred from 
the displacemeat required at the supporting beam to return the element to 
null. In the present configuration, the supporting beam is driven by a 
differential micrometer powered by a small variable-speed motor. Two 
independent measurements were made of the beam motion. The counter shown 
in Fig. ll measured the rotation of the micrometer drive shaft and counted 
in units of 10°76 inches of axial displacement. In addition, a Schaevitz coil 
was mechanically linked to the beam through the drive wire. The null position 
of the element was monitored by a second Schaevitz coil, as in the original 
design. The demodulated output from the Schaevitz coils was low-pass filtered 
with a time constant of 0.25 sec. The dashrots shown in Fig. 11 were filled 
with Dow Corning 710 silicone oil having a viscosity v = 5.0 em? /séc. A 
thermocouple measured the temperature of the balance chamber. 

The rectangular floating element, shown in Fig. 12, fs 0.622 cm in 
the streamise direction and 3.785 cm in the cross-stream direction. The area 
of the element is thus A = 2.356 cm, The gap is 0.007 cm upstream and down- 
‘stream when the element is nulled, and 0.010 cm on erch side. The element was 


flush with the surrounding surface within 0.G00l cm. No correction was made 


to account for the effect of the gap on the measured force. 
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The balance was calibrated using the technique described by 
Coles (1953). The beam displacement required to rerurn the element to 
null was measured with the balance tilteu at various ang'es with respect 
to the horizontal. These measurements were repeated after adding varior. 
small weights to the element, and the results were analyzed to yield the 

ass of the umweighted element and the spring constant for the ‘lexures. 

Four angles were used between O66" and 0.6, with cottnee of 0, 5, 10 and 

20 grams. The spring constant was measured to be 73.98 gm/cm, with a maximum 
deviation of 2 parts in 1000. 

The JPL 20-inch tunnel is a variable-density facility. The balance 
was located on the tunnel ceiling, which flexes with changes in iree-stream 
static pressure. To compensate for the resulting zero offset in the surface-_ 
friction balance, the element was covered by a thin gasket-sealed plate, which 
was held in place by evacuating the balance chamber, and flow was established 
at the desired Mach number and Reynolds number. When the balance achieved 
thermal equilibr the null position of the element was measured. The fle 
was then bypassed, e element cover was removed, and flow was re-established 
at the same free-s ~u conditions. When thermal equilibrium was again achieved, 
the null position of the element was again measured. The difference between 
the beam positions required for null, with and without applied shearing stre“s, 
is a direct measure of the force on the element. 

The measurement just described is not quite correct, because any 
differential pressure between the sealed balance chamber and the test section 


causes additional bending of the b:lance port. The resulting zero offset was 


determined with the tunnel off ard with the test-section pressure set at appre- 


priate values. The maximum correction applied to the data wis 6.5 percent. 
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c® 
Finally, beccuse the equilitrium temperature may not be the same for the 


various null measurements, effects of thermal distortion must also be con- 
sidered. The zero offset from this source was measured separately; the 
maximum correction applied to the data was 0.4 percent. The streamvise 
force on the element due to free-stream pressure gradient is negligible. 
Measurements using the balance were made at only one station, 
as indicated in Table 1. The data, together with the free-stream static 
and stagnation contitions, were recorded by the data-acquisition system, 
The surface-friction balance yields a direct measurement of the tangential 


stress on the floating element, 


F 
T ran (50) 


are listed in Table 3. 

B. Preston Tube 

An independent estimate of surface friction was obtained from 
measurements with a Preston tube, a flat-faced circular cylinder in contact 
with the wall. For the JPL experiments, three probes were used, having 
outer diameters D of 0.082 cm, 0.162 cm, and 0.317 cm, with a ratio of inner 
to outer diameter of 0.60. The largest probe was tested only at station 


JPL-2. The probes were positioned at the wall by the traverse mechanism. 


We 
Waar oo 
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For the CIT experiments, a single Preston tube wa: used, having an outer 
diameter of 0.210 cm, with a ratio of inner to outcr diameter of 0.76. 
There is no concensus on the question of proper interpretation of 
Preston-tube data at supersoric speeds. For example, Hopkins and Keener 
(1966) took as the geometric parameter 
0,4,D 


Re, = , 52 
eee (52) 


and as the response parameter, the Mach number a implied by the ratio of 
Preston-tube pressure to local static pre.sure. They proposed, as a cor- 
relation between this parameter HM, and the surface friction coefficient 


Ce» the expression 


uw \? 1.132 
£.(T’) Re7( 2} = 32.885 iP (T’) Re” c | (53) 
2 p \H, : u‘2 D “Et , 
where 
> 2 ’ 
’y a [ & ° 
£,(T') (“s) ae (54) 


with »’ and 9’ evaluated at the reference temperature T’ proposed by Sommer 


and Short (1955), 


Ee 140.035 4,2 + 0.05(s* - 1) (55) 
Tt, 7 e T, 7 


Values of local friction coefficient Ce calculated from these 


equations are included in Table 3 and are listed in more deta!1l in Table A2 


of the Appendix. 
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Bradshaw and Unsworth (1974) have taken the position that only 
wall quantities should appear in any Preston-tube correlation and that it 


is unrealistic to insist on an explicic formla for C They took as a 


f° 
point of departure a recent survey by Allen (1973). Using Allen's own 
calibration data (but not other data considered by Allen), they proposed 
a formula representing Patel's low-speed calibration (1965), with an addi- 


tive term to account for compressibility. The forma applies for adiabatic 


flow and for 50 <D* < 1000, where 


, (56) 


and Te Pai as before. 

Unfortunately, Allen has recently reported (1977) that his published 
friction measurements are incorrect, because of a defective or poorly cali- 
brated balance. He did not repeat h’ experiment, but simply replaced the 
original measured values of surface friction by computed ones. He also 
made the corresponding revision in the Bradshaw-Unsworth formula (again 


using only his own revised calibration), with the final result 


- s at : 
ae 96 + 60 tee, (35) + 23.7 10810 (35) 


0.30 
+ 10° m7 [(o*) - 2.20| (57) 
where 
fs 2 
2 T 

um? 2 : (58) 

T vRT,, 

30 


AEOC-TR-78-21 


and 
ci = 2 zt (59) 


Here P, is the pressure measured by the Preston tube and P, is the ambient 
static pressure. 

The present Preston-tube data have also been processed in terms 
of Eq. (57) to obtain the values of local friction coefficient C, which are 


£ 
listed in Table 3 and in Table A2 of the Appendix”. 


C. Friction-Law Scaling 

It has been pointed out by Spalding and Chi (1964) and others 
that most analytical formulations for compressible turbulent boundary layers 
are reducible to a description of the surface friction in terms of the surface 
friction for an equivalent incompressible boundary layer at a differeni 
Reynolds number. The equivalence is usually expressed by two semi-empirical 


and F_, thus: 


scaling functions Fe 8 


C, = F.C. , (60) 


and 


i 
Reg F Reg . (61) 


*the experiments described in the present report should eventually be viewed 
as a producer rather than as a consumer of Preston-tube calibration data. 
However, a full-scale critique of the Preston-tube technique for supersonic 
flow is outside the scope of this research, and the present measurements 


have therefore been interpreted as if no friction-balance data were available. 
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For the particular case of adiabatic flow, che scaling functions 
= 2 
Ee See a ep ae haa [=F , (62) 
3 -1 T,\t/2]? (- n’) sin a 
sin [1 - Tr 


and 


, (63) 


@D 
& |a™ 


were first proposed by Wilson (1949). Use of Eq. (63) for flow with heat 
transfer was later recommended by Van Driest (1955) and is sometimes re- 
ferred to as Van Driest I1. For adiabatic flow, given M, and r (and hence 
Fe and Fy)» and given also Reg. 
puting Re a" from Eq. (61) by looking up the associated value for ec us ing 


the surface friction is determined by com- 


some convenient low-speed friction law; and finally by computing ce from 
Eq. (60). That is, 
1 4 i 1A 
Cc, Fe Ce (re; => Ce (rg Reg] . (64) 
A “convenient low-speed friction law" is implicit in the survey 
by Coles (1968), who recommends the value Ii = 9.62 for flow at constant 


> 5000. For lower Reynolds mimbers, we multiply 


pressure when Re 
M(6u,/v) from Table 2 of Coles (1962, .ppes:ix A) by 0.62/0.55,. With 

i _ 2 
bu_/v as independent variable, the quantities uo/e, and Ce =2 (u,/a,) 
follow immediately from the local friction law (21) above. Finally, we 


compute Re gs and Ret from Eqs. (7) and (8) of Coles (1968), after replacing 


8 


Re gx by Re gx - 65 to take account of the real profile in the sublayer. The 
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result of these calculations is recorded in Table 4. 


Values calculated for C, from Eq. (64), using futerpolation tn 


f 
Table 4 to define the function c*(Re gy). are included in Table 3. 
V. Discussion and Conclusions 

According to Table 3, five different methods have been used to 


measure or to estimate the local friction coefficient C, for the present 


f 
experiments. Three of these methods (Van Driest scaling, Preston tube, 
friction-law scaling) depend on some empirical means for taking account of 


compressibility. In particular, the friction-law scaling of Section IV-C 


e by Res and the measured or estimated Ce by Gr 


When the data of Table 3 are subjected to this same scaling, they appear 


replaces the measured Re 


as shown in Fig. 13. The solid curves represent low-speed experience 
according to Table 4, Any discrepancy between the data and the solid curves 
should not be interpreted as error, because the friction-law scaling itself 
would then have to be viewed as error-free. This scaling is in fact of 
uncertain accuracy, and is used primarily to remove most of the effects 


of Mach number in the data, so that one technique for evaluating C, can 


f 
be readily comparei with another. The required displacements from (Cy. 
Re ,) to ce Re) are indicated by the Line segments rext to the lowest 
curve in Fig. 13. These displacements are not very substantial, because 
the Mach numbers for the present experiments are relatively low. 

We consider the floating-element friction data to be the most 
reliable data in the figure. For the estimates of Ce from dA/dx, the 


scatter is large, as expected. The largest scatter, however, is in the 


Preston-tube data, indicating that this technique needs further development. 
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Of the two Preston-tube correlations, the one by Bradshaw and Unsworth 
(as revised by Allen; see Section IV B) underestimates Ce slightly, 
especially at M, = 2.2. The correlation by Hopkins and Keener is satis- 
factory except at M, = 2.2, where it overestimates Ce by a large amount. 
The profile fic gives values for Ce which are systematically a little high 
(except for the CIT measurements), with a maximum discrepancy of about 
6 percent at M, = 2.2. 

In general, scaling of the measured mean-velocity profiles 
according to the Van Driest version of the mixing-length theory (Eq. (18)) 
scems to be quite successful. The scaled profiles can be well represented 
by conventional low-speed wail-wake similarity laws. Except for discon- 
tinuities in slope at the edge of the boundary layer, the inferred profiles 
for v/u and Th provide a quite acceptable standard for interpretation 
of the LDV measurements reported in Part II of this work. Estimates for 
v/u and tit based on the assumption of similarity in y/& or y/8 are less . 


satisfactory. 


34 


te ey ee 


AEOC-TR-78-21 


References 


Abbiss, J. B. 1976 Development of pheton correlation anemometry 


for application to supersonic tlow. AGARD Conf. Proc. CP 193, 


“Applications of Non-intrusive Instrumentation in Fluid Flow 
Research," Paper il. 

Allen, J. M. 1973 Evaluation of compressible-flow Preston tube 
calibrations. NASA TN D-7190; see also AIAA J. 11, 1461-1462, 
1975. 

Alleu, J. M. 1977 Re-evaluation of compressible-flow Preston tube 
calibrations. NASA T X-3488. 


Bradshaw, P. and Unsworth, K. 1974 Comment on "Evaluation of Preston 


—— ee 


tube calibration equations in supersonic flow," with Reply by 
J.M. Allen. AIAA J. 12, 1293-1296. 


Bushnell, D. M. and Morris, D. J. 1971 Shear-~stress, eddy-viscosity, 


and mixing-length distributions in hypersunic turbuleat_ boundary 
layers. NASA TM X-2310. 


Coles, D. 1953 Measurements in the boundary laver on a smooth flat 
plate in supersonic flow, Instrumentation and experimental 
techniques at the Jet Propulsion Laboratory. JPL Rep. No. 20-70. 

Coles, D. 1962 The turbulent boundary layer in a compressible fluid. 
Rand Corp. Rep. R-403-PR. 

Coles, D. 1968 The young person's guide to the data. Proc. 1968 AFOSR- 
IFP-Stanford Conf., "Computation of Turbulent Boundary Layers," 


Vol. Ii (D. Coles and E. Hirst, eds.), 1-45. 


35 


AEOC.TR.78-21 


9. Danberg, J. FE. 1971 A_re-evaluation of zero pressure gradient 
compressible turbulent boundary layer measurements. AGARD 
Conf. Proc. CP-93, “Turbulent Shear Flows," Paper 1; also 
BRL Rep. No. 1642, 1973. 


10. Fenter, F. W. and Stalmach, C. J., Jr. 1957 The measurement of 


local turbulent skin friction at supersonic speeds by means 

of surface impact pressure probes. Univ. Texas, Austin, 

Rep. DRL-392 (CM-878); see also J. Aero/Space Sci. 25, 793-794, 
1958. 

ll. Fernhu z, H. H. 1976 Compressible turbulent boundary layers. In 
VKI Lecture Series 86, "Compressible Turbulent Boundary Layers," 
Vol. I. 

12. Fernholz, H. H. and Finley, P. J. 1977 A critical compilation of 
compressible turbulent boundary layer data. AGARDograph No. 

223. 

13. Horstman, C. C. and Owen, F. K. 1972 Turbulent properties of a 
compressible boundary layer. AIAA J. 10, 1418-1424. 

14. Hopkins, E. J. and Keener, E. R. 1966 Study of surface Pitots for 
measuring turbulent skin friction at supersonic Mach numbers-- 
adiabatic wall. NASA TN D-3478. 

15. Johnson, D. A. and Rose, W. C. 1973 Measurement of turbulence transport 
properties in a supersonic boundary-layer flow using laser velocimeter 
and hot-wire anemometer techniques. AIAA Paper 73-1045; see also 
AIAA J. 13, 512-515, 1975. 

16. Klebanoff, P. S. 1954 Characteristics of turbulence in a boundary 
layer with zero pressure gradient. NACA TN 3178; also TR 1247, 1955. 


36 


“4 
‘ 


17. 


18, 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


AEDC-TR-78-21 


Kleinstein, G. 1967 Generalized law of the wall and eddy-viscosity 


model for wall boundary layers. AIAA J. 5, 1402-1407. 


Liepm-nn, H. W. and Ashkenas, H. 1947 Shovk-wave oscillations in wind 


tunnels. J. Aeron. Sci. 14, 295-302. 


Maise, G. and McDonald, H. 1968 Mixing length and kinematic eddy 


viscosity in a compressible turbylent boundary layer. AIAA J. 6, 
73-80. 


Meier, H. U. and Rotta, J. C. 1970 Experimental and theoretical 
investigations of temperature distributions in supersonic 
boundary layers. AIAA Paper 70-744; also AIAA J. 11, 2149-2156, 
1971. 

Michel, R., Quemart, C., and Elena, M. 1969 Distributions de vitesses 
des couches limites turbulentes en écoulement compressible, 
uniforme ou accélérg. La Recherche Aérospatiale, No. 128, 33-47. 

Moore, D. R. and Harkness, J. 1964 Experimental investigation of the 


compressible turbulent boundary layer at very high Reynolds numbers, 
AIAA Paper 64-592; also AIAA J. 3, 631-638, 1965. 


Morkovin, M. V. 1961 Effects of compressibility on turbulent flows. 
Proc. Colloq. "Mecanique de la turbulence," CNRS (1962), 367-380 
(Proc. reprinted as "Mechanics of Turbulence," Gordon and Breach, 
1964). 

Patel, V. C. 1965 Calibtration- of the Preston tube and limitations on 
its use in pressure gradients. J. Fluid Mech. 23, 185-208. 


Rotta, J. C. 1960 Turbulent boundary layers with heat transfer in 


compressible flow. AGARD Rep. 281. 


ORIGINAL BAGE 8 
OF POOR QUALITY 


37 


AEDC-TR-78-21 


26. Sandborn, V. A. 1974 A_review of turbulence measurements in 
compressible flow. NASA TM X-62337. 

27. Sommer, S.C. and Short, %. J. 1955 Free-flight measurements of 
turbulent boundary-layer skin fri fon_in the presence of 
severe aerodynamic heating at Mach numbers from 2.8 to 7.0. 


NACA TN 3391; see also J. Aeron. Sci. 23, 536-542, 1956. 


28. Spalding, D. B. 1961 A single formula for the “law of the wall". 
Trans. ASME 28E (J. Appl. Mech.), 455-457. 

29. Spalding, D. B. and Chi, S. W. 1964 The drag of a compressible 
turbulent boundary layer on a smooth flat plate with and 


without heat transfer. J. Fluid Mech. 18, 117-143. 


30. Squire, L. C. 1971 Eddy viscosity distributions in compressible 
turbulent boundary “ayers with injection. Aeron. Quart. 22, 
169-182. 

31. Sturek, W. B. 1973 Calculations of turbulent shear stress in 
supersonic turbulent boundiry layer zero and adverse pressure 
gradient flow. AIAA Paper 73-166. 


32. Van Driest, E. R. 1951 Turbulent boundary layer in compressible fluids. 
J. Aeron. Sci. 18, 145-160, 216. 


33. Van Driest, E. R. 1955 The turbulent boundary layer with variable 
Prandtl number. "50 Jahre Grenzschichtforschung" (H. Gortler 
and W. Tollmien, eds.), Vieweg, Braunschweig, 257-271; see also 


The problem of aerodynamic heating, Aeron. Eng. Rev., Oct. 1956, 
26-41, 


38 


AEOC-TR-78-21 


34. Wieghardt, K. 1943 Uber die Wandschubspannung in turbulenten 

Reibungsschichten bei verdnderlicnem Aussendruck. KWI, 
Gottingen. U & M 6405; see also K. Wiegharde and W. Tillmann, 
“Zur turbulenten Reibungsschicht bei Druckanstieg," KW1, 


Gortingen, U & M 6617, 1944 (translated as “On the turbulent 


friction layer for rising pressure,'’ NACA TM 1314, 1951). 


35. Wilson, R. E. 1949 Turbulent boundary-layer characteristics at 


supersonic speeds ~ theory and experiment. Univ. Texas, Austin, 
Rep. DRL-221 (CM-569); see also J. Aeron. Sci. 17, 585-594, 1950. 


36. Winter, K. G. and Gaudet, L. 1970 Turbulent boundary-layer studies 
at high Reynolds numbers at Mach numbers between 0.2 and 2.8. 
ARC R & M 3712 (1973). 


37, Yanta, W. J. and Lee, R. E. 1974 Determination of turbulence transport 


properties with the laser Doppler velocimeter and conventional 


time-averag-d mcan flow measurements at Mach 3. ALAA Paper 74-575; 
see also AIAA J. 14, 725-729, 1976. 


FI 
ORIGINAL PAGE 
OE POOR QUALITY 


AEDC-TR-78-21 


Station 


Table 1. 


Station 


Table 2. 


Location 
(em from floating- 
element balance) 


ny 
vey 
3 
H 
a 
° 
u 
tt 
Q. 


High-Speed Flow Measurements 
JPL 20-Inch Wind sunnel 


Location 
(cm from leading 
edge) 


Pitot Tube 


4 
9 
4 
4 
.6 
8 
el 
3 
-6 


Low-Speed Flow Measurements 
CIT Merrill Wind Tunnel 


40 


Preston Tube 


Pitot Rake 


Balance 


Preston Tube 


xx MK KR KK KK 


AEDC-TR-78-21 


Se 6luVU* eS6Tuu* 

£56100" 276100° vabTOO" Oc oT uu" ¥vLoTOUL® 276100 °* 
€S6100* L86100° 

7L6100° SE6TO0* T26T 00" £66T00° 

H65TOL® s00Z00° 

€0TZou" SOT Zou" 

9TTcuV® $3020 9° £2¢200° 9902 GU" Ocl2é0ou* 270200 * 
LTT¢< 00° cETZOU® 

ByTc00° V6L2 00° 6€ Te uu" 6UT 200° 

Litéuv* 9ETCOU® 

020200" 1970200° 

2€0200° %766100° S86T00°" €¥ol00® $90Z00° 266TU0* 
S2ucuu* 9S0200° 

80200" sT¢¢00° 2TbZ@Uu" £50200" 

eL0cOu® 06020uU° 

261c 00° S61T200° 

O6Tc0U? s9TZ00° 901z200° 60T2 60° BoT 200° 96020U° 
%61T200° 961200° 

L£22200° 691200° €Z12U0° T02200° 

6422 00° 6L1200° 

£89¢00° $45 200° £07200" GL9C00° 

Se@LZ00° vOLZ 00° TT9Z00° 689200° 

6¢€2200° $28200° 9ELZUO® 689200° 

eeulzZ0u® . ZLLZ00° ECBIZLO® 954200° 89L200° 
bO8zco00° 2£9Z200° €9S200° 982L7200° 

92820U° ST2L200° O€ ¥20U" £824 200° 

T2€e00° 286200° L8UE00® 

£0L€00° tTece00* 8S7e00* 

01 700° Tt6¢00d° 4$0%700° 

(n/@) QW/H) 


QsiNgwO)  32NV1V98 NOLS3ta NOLS3ud XO10 #2 


a2 


ABVWWHS NOITLIIYS NIWS 


*€ JIGVL 


"OQUu9e4y 
“0001s 
*U6TC 
“O9ELE 
*O96EE 


*OQLSve 
*Olzte 
"U7GE¢ 
*Osblc 
*uLlol 


*OcoLe 
“OLYSE 
*OBCLE 
*OCE*E 
"Ovrle 


“QUE? 
*00%22 
*0oT22 
°08 102 
“OL38T 


*8908 
°SLoL 
*OLeL 
°%099 
*6029 
*2t6S 
*86L2 
*slst 
*6zOT 


VL3HI-34 


oloz® 
lcol°* 
Ov76L° 
€76L° 
OfoL*® 


SoboL*® 
STU? 
o7Us® 
cuyl® 
8Sol°® 


Stos°* 
T€6S ° 
¢S6S° 
2496S ° 
tLoS° 


296S° 
8glOv® 
986s° 
£cos° 
Lets? 


oLot® 
csol* 
9EOt? 
TEOT® 
cLor® 
BsuoT® 
osot® 
osuolt® 
Qsaol? 


3W 


S-Vdf 
9-1df 
€-1d 
etd 
T-16f 


S-i\df 
9-Idf 
t-1df 
é-1df 
T-1df 


s-df 
9=-10F 
e-1df 
¢-14r 
t-1d6¢ 


s-T1dt 
9-100 
€-1d 
e-idt 
i-l0r 


6-119 
G-119 
z-:1) 
9-119 
S-1ID 
teLID 
€-L119 
en-L13 
t-119 


NOILVLS 


BAGE IS 
ALITY 


ORIGINAL 
OF POOR QU 


41 


| 
| 
| 


L74100° 
95%100° 
297100° 
9L%( 00" 


€¥S10vu° 
O6Stou* 
€O0vtTOo* 
L09TO0O* 


2e21uv® 
S6LlO0vu* 
SeL1u00* 
2081 Ou" 


leo tou" 
¥S6TO0O° 
#S6L00°" 
766100° 


764 100° 
$061 00° 
$.6T00° 
0€6100° 
£96100° 


290¢00° 
920¢00° 
9240200" 
£OTc0u* 
9971TZ00° 


AEDOC-TR-78-21 


OQ3LNdwWOdD 


$74100° 


ees too° 


BuLTOO° 


2498T00° 


£751 00° 


2502 00° 


4 


Ble Too® 


SBE T UO" 


Lo%Tu0* 


eLe100* 


£691 00° 


102 T 00° 


gOwToO* 


£6LT 00° 


OT8TO0° 


€9¥ 100° 


¥00200° 


720200° 


(t/@) 
NOL$3ud 


32 


€zgt00° 


€tgtoue 


CvITOC® 


OvLT 00° 


20xT0V0" 


BLLTOO" 


€ToTUU® 


906T OU? 


OLvTUu® 


2eotuvu® 


Twuz00° 


811200° 


(3/H) 
NULS 34d 


(°iN09) 


405100° 
2estou° 
O€S TOU" 
7S T0U% 


929100° 
€t 9100" 
699 TOO" 
9s9ICO° 


ZtB TOU" 
ug¥TOuU® 
¥SBTO0° 
448 TOU" 


6$6100° 
€¥6t00* 
€v6Tou® 
0v0200* 


TL6TOL* 
S¢6100° 
€$6 100° 
076 100% 
O26 T00* 


490200° 
T80200* 
£60200" 
$9020u* 
$01z00° 


ilj 


*€ +evl 


99% 100° 


Z2€ST00° 


0Ss2 100° 


9S8TOO* 


710200° 


%$0200° 


"O90EY 
"O09 T? 
°0LS0% 
*Osuse 


*090S2 
*069%2 
*O2se2 
*OLOE?S 


"Ol¢0 
*QU6LE 
"OSSLE 
*OECLE 


*OoT%2 
*O06t? 
*OseTZ 
*ObL6T 


"OSSI% 
*OO66E 
*OOSBE 
*OSZ9E 
“OCECE 


*OSBE?2 
“O98 22 
*O2L2e 
* 06802 
*0s98t 


VL3SHL-3u 


Loit??é 
ocet *< 
LELT°? 
clgt*2 


fcLetc?e 
27918? 
999TPS 
C220 °? 


veTtet 
sclect 
ezle*l 
ceuE Pl 


Tstert 
Lote rt 
SIcect 
Tote °T 


9096" 
Le96° 
E196° 
9¢296° 
8946" 


1s96° 
2L96° 
6TL6° 
6996° 
97976° 


aw 


s-1d0f 
7-f 
€-1df 
g-106f 


$-1¢r 
9-100 
€-0P 
e-0f 


S-16f 
7-Vdl 
€-1dF 
é Wf 


GS-16 
9-10 
€-106F 
co idt 


$-100 
9-10F 
€-IdF 
é-10f 
t-1df 


$-T4r 
9-160 
e-1df 
2-160 
i-U0r 


NOILVIS 


AEDC-TR-78-21 


Table 4. 


tAw-SPFED FRICTION LAW 


NELTA-PLUS PT WESUT ce RE-NSFARQ RE-THETA 
24M. ennn 18,37 e9NSOPR ASN. 4340, 
300. 2138 19.57 20052771 R99. anT. 
4nn, 2759 20.88 2ON45AR 1794. R96. 
500. 0 438 F1L-RI 2 INS204 1497. 11989. 
60, 240A 272 5A e9NIAQA77 2127. 149A, 
RN0, 2 GAS 23-67 2 J03570 2AH2Z «© 7114. 

19000, 294) 24249 2003335 3R24, 2749. 
159n, 2598 25275 eON3N1H 5°10. LANA, 
2000. 2479 Phe 56 e007 R36 TAT. BRAS, 
aonn. ean 272558 2997 63486 11920, RANK, 
4nnn,. 2679 PR?5 enn7 505 15270. 11980, 
§50n0, «6An 2R,R0 20997417 19820. 159469. 
4000. 2fh2n 29.2% 09097339 23770. TRIED, 
Rann, 2629 279,94 2907730 ZLARO. P4IAN, 
19000. 262N 302.49 29971527 395R0. 3NKGN, 
15000. 2670 31.48 2007018 59330. GA3A0. 
20000. e670 32.18 2001931 79090. 62210. 
30000. 2A? 33.17 eNAYRIA 11AA00. INTO, 

_ pack & 

oricit a Qu aT 

oF PO 
43 


AEOC-TR 78-21 


2.4 


° 
e 
3 
8 
Me=1,32 
-—-—_ —-—_ + —--- —e--- B- S22 $4 4 os ti, > _——__—__ 
° 
e 
° 
e 
MNe= .97 3.0 : 
ooo ~~ eg ~--- * of O80 — A 
t 
Me= .80 mes 
Oe ee ee >---B- gee eo — BR. Tage a SE ee 
e 
° 
te 
Me= .60 
_—_———— Sa en egg et a ef 
* ° 
ae 
> 
t=) 
— 
[ Cnn COD Cans Sey COED Sens CS SaGeeey Games Wn Wan ry ra : a oe a ae 


-1NO -3CG -100 ~-H9 -60 -4O) -@u 0 eg ag 60 BO 


X cm 


Finure 1, Free-stream Mach-Number Distribution, Re, =40,000, 
JPL 20-inch Wind Tunnel 


44 


- 


AEOC-TR-78-21 


150 


0125 


M Mes .11 


—9-8—_9—_9 9 ——_-g§——_- _——__ 


2100 
2075 
.0SO } 


025 


LOV 
- TE 


«000 
so 100 150 200 250 


X (cm) 


Figure 2. Free-stream Mach Number Distribution, 
CIT Merrill Wind Tunnel 


B lb 
IGINAL PAG 
OF POOR QUALITY 


4§ 


: e 


AEOC-TR-76-21 


Q 


= 
{ 
Q 
~” 
R Q 
Q 
toe, Z mt 
coe, “roel. Qa 
~ 
Q 
~_ 
o 
ioe} 
~ 
wa 
ae wn 
od 
mM 


o a a@ 84 a4 


es cee ee + Sao ee, See eee es 


“WLSHL/A 


U/UE 


Figure 3. Mean Velocity Profiles. 


Re, = 23,000 


46 


AEOC.TR-78.21 


| | 
Poe hts ears (ea 
| 


1 


Mos el 
"Rea 
get 


d 


Q 
~ 
Je 
~ 
Jo 
~” 
3 
2 ras 
iui o 
¢ 
es r 
~ o 
fo) fod 
~ 
o 


6 


6 
U/UE 


Figure 4. Mean Velocity Profiles. 
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Figure 11. Floating-Element Balances. 
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Figure 12. Floaving-Element Balance. 
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Appendix 
Tabulation of Experimental Data 


Table Al 


AEOC-TR-78-21 


Table Al sunmarizes various parameters for the profiles obtained 


from Pitot-tube surveys at the stations listed in Tables 1 and 2 of the 


main text. The quantities uss Tl, and 6 are derived from the least-squares 


fit of each profile to the wali-wake formla, 


* yu 
Be bin tee +2 Fin’ (FX) : 
uo wv x 26 


where u™ is velocity scaled according to Van Driest, 


u* -{(2y” du. 


The two integral thicknesses 6* and @ are defined by 


6 
* . 
ie 
and by 6 
ae, Lae (: - 2) ay . 
o Pere Ve 
The quantity H is che ratio 6*/0. 
PAGE bb 
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(11) 


(4) 


(5) 
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Table A2 


Tabl= A2 contains Preston-tube data, including probe diameter D, 


pressure coefficient 


p 2 7? (59) 


vP,Me 


ai friction coefficient, as inferred by two different correlation methods. 
The method of Hopkins and Keener (1966), denoted by H/K, uses the formla 
mM \2 1.132 
£,(T’) Rea (—2] = 32.885 [« (T’) Rez c | 
2 p\™ : 2 pd “¢ , (53) 
which is explicitly soluble for Ce. The quantity 4, is the Mach number 
implied by co and the quantities Re) and £,(T’) are defined by Eqs. (52) 


and (54) of the main text. The method of Bradshaw and Unsworth (1974) as 


revised by Allen (1977), denoted by B/U, uses the formla 


: f) 20 fina) 
c = 96 + 60 login (3)+ 23.7 [sexo 50 


(57) 
4 +,0.30 
+ 10° lo yeren e 2.30] ; 


2 


This formula is soluble for Ce only by iteration, since the quantities pt 
and M_ both depend on Cy. These quantities are defined by Eqs. (56) and 


(58) of the main text. 
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Table A3 


Table A3 includes estimates of the pressure-gradient parameter 


1 4P 
BETA eng 72, (Al) 


DTDX = de ‘ (A2) 
These are assigned to stations CIT-6 or JPL-4 as appropriate. The quantity 


, (A3) 


where Ce is the value measured using the floating-element balance, should 


equal unity if the experiment is free of error. 


Tables A4-A14 


Tables A4-Al4 and the associated Figs. Al-A44 are a detailed re ord 
of data obtained from the Pitot-tube surveys. Each table heading incluves 
the integral properties 6* and @ for the profile, together with the friction 
coefficient Ce measurec using the floating-element balance, where applicable 


(1.e., station JPL-4). Listed neat are the profile parameters u_, il, 6, from 


the profile fit, wa... the associated friction coefficient 


2 
L (=) 
Cc. 2 De m" * (22) 
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The range of y used in the fit is specified as YMIN, YMAX. The mean square 
deviation of the fitted data from the wall-wake formula is given as CHISQR. 
The variable for this calculation is the Van Dricst velocity u*/u,. 

The body of the tables lists the distance from the wall as y, as y/6, 


and as 


Y-PLUS =——_.,, (A4) 


Also listed are the local Mach number M, density p, and velocity u (all 
normalized by the corresponding free-stream values), and the Van Driest 
velocity 


* 
U-PLUS =— . (A5) 


Finally, the tatles give the shearing stress By Le and the normal velocity 


v/u computed from 


an-4p 

~~ 21 ve 37 
T “\ 2Q9 - P : (37) 

w e e 

and 
T ou 

v w ee P 

Fs a oa, (38) 
4a Pee pu (2Q, P.) 


where P and Q are definite integrals defined by Eqs. (44) and (45) of the 


main text. 
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Figure A1, Mean Velocity Profiles. 
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Figure A4, Shear Stress Distribution. 
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Figure f&-S. Mean Velocity Profiles. 
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Figure A41, Mean Velocity Profiles. 
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Figure A44. Shear Stress Distribution. 
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Figure A41, Mean Velocity °Srofiles,. 
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Figure A44, Shear Stress Distribution. 
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Nomenclature 
Symbol Equation Meaning 
c (19) constant in wall law (5.0) 
Ce (7) local friction coefficient 
q (59) pressure coeffictent for Preston 
tube 
D (52) diameter of Preston tube 
f(y") (40) function in wall law 
£,(T’) (54) function of reference temperature 
Fe (62) scaling function for Ce 
Py (63) scaling function for Rey 
u (6) boundary-layer profile form parameter 
n (14) function of Mach number 
M, (53) pressure Mach number for Preston tube 
M (58) friction Mach number for Preston tube 
MOMB (A3) measure of momentum balance 
P,Q (30) ,¢€31) definite integrals of velocity profile 
r (2) temperature recovery factor (0.885) 
Rey (22) Reynolds number Lased on D 
Beg Reynolds number based on § 
u,v streamwise and normal velocity 
components 
u,, (15) friction velocity 
U (25) dimensionless velocity scaled according 


to Van Driest 
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Subscripts 


O, 
©), 
or 


Equation 


(20) 


(Al) 


(4) 


(5) 
(9) 


(19) 


Superscripts 
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eh 
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Nomenclature 


+e bd 
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Cont. 


Meaning 


streamwise and normal coordinates 
distance from wall in outer variables 
pressure-gradient parameter 
boundary-layer thickness 
boundary-layer displacement thickness 
boundary-layer momentum thickness 
Karman constant (0.41) 

viscosity 

kinematic viscosity 

strength of wake component 


shearing stress 


edge or external value 
stagnation value 


wall value 


value for incompressible flow 
value at effective temperature 


value made dimensionless with 


ul, v 
T OW 


value for Van Driest scaling 


257 


